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Abstract 



In recent years, there has been considerable interest in theories for- 
mulated in anti-de Sitter (AdS) spacetime. However, AdS spacetime fails 
to be globally hyperbolic, so a classical field satisfying a hyperbolic wave 
equation on AdS spacetime need not have a well defined dynamics. Nev- 
ertheless, AdS spacetime is static, so the possible rules of dynamics for 
a field satisfying a linear wave equation are constrained by our previous 
general analysis — given in paper II — where it was shown that the possible 
choices of dynamics correspond to choices of positive, self-adjoint exten- 
sions of a certain differential operator, A. In the present paper, we reduce 
the analysis of electromagnetic, and gravitational perturbations in AdS 
spacetime to scalar wave equations. We then apply our general results 
to analyse the possible dynamics of scalar, electromagnetic, and gravita- 
tional perturbations in AdS spacetime. In AdS spacetime, the freedom (if 
any) in choosing self-adjoint extensions of A corresponds to the freedom (if 
any) in choosing suitable boundary conditions at infinity, so our analysis 
determines all of the possible boundary conditions that can be imposed 
at infinity. In particular, we show that other boundary conditions besides 
the Dirichlet and Neumann conditions may be possible, depending on the 
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value of the effective mass for scalar field perturbations, and depending on 
the number of spacetime dimensions and type of mode for electromagnetic 
and gravitational perturbations. 

1 Introduction 

Anti-de Sitter (AdS) spacetime 1 has come to play a central role in a number of 
contexts in theoretical physics [1], most prominently in the "AdS-CFT" corre- 
spondence [2, 3, 4] (see also e.g., [5] and references therein). Although AdS space- 
time is maximally symmetric and geodesically complete, it admits no Cauchy 
surface and thus provides a very simple but nontrivial example of a non-globally- 
hyperbolic spacetime. Consequently, it is far from obvious, a priori, that one can 
define a sensible, deterministic field dynamics in AdS spacetime, even for classi- 
cal, linear fields: If a classical field satisfies a linear hyperbolic equation, then it is 
locally uniquely determined by suitable initial data on a spacelike hypersurface; 
however, in order to guarantee that a solution corresponding to given initial data 
exists globally and is globally unique, it is essential that the spacetime be globally 
hyperbolic and that the data be specified on a Cauchy surface (see, e.g. [6]). 

For the case of a scalar wave equation in AdS spacetime, it is clear that— 
at least for some values of mass and coupling to curvature — it is necessary to 
specify boundary conditions at infinity in order to obtain a well posed, deter- 
ministic evolution law. This issue has been investigated by a number of authors 
(see, in particular, [7, 8, 9]) and some proposals — such as Dirichlet or Neumann 
conditions — have been made. However, it is far from obvious that the boundary 
conditions that have been proposed encompass all the possibilities for dynamical 
evolution laws in AdS spacetime. New dynamical evolution laws would provide 
new "bulk theories" in AdS spacetime, and could be of interest in view of the 
AdS-CFT correspondence 2 , among other reasons. 

The purpose of this paper is to apply our general analysis of dynamics in static, 
non-globally-hyperbolic spacetimes [16] to systematically determine all possible 
boundary conditions that can be imposed at infinity for scalar, electromagnetic, 
and gravitational perturbations of AdS spacetime. In our previous paper [16], 
we proved that any dynamical evolution law for a scalar field in a static, stably 
causal spacetime that (i) locally agrees with the wave equation, (ii) admits a 
suitable conserved positive energy, (iii) is compatible with the time translation 

x In this paper, by m-dimensional AdS spacetime of curvature radius £, we mean the universal 
covering space of the hyperboloid composed of points whose squared distance from the origin 
is — £ 2 in R m+1 with a flat metric of signature K..+. 

2 Some implications of different choices of boundary conditions for AdS-CFT correspondence 
have been discussed in [fO, 11]. In particular, it was pointed out that multi-trace deformation 
of a boundary quantum field theory can be incorporated in the AdS-CFT correspondence by 
using a general boundary condition for the corresponding bulk field [12, 13] (see also [14, 15] 
and references therein). 
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and time reflection symmetries, and (iv) satisfies a certain convergence condition 
must correspond to a prescription of the type first proposed in [17]. Consequently, 
the general prescription for a dynamical evolution law for a scalar field in a static, 
stably causal spacetime arises as follows: 

Let (M, g^) be a stably causal (but not necessarily globally hyperbolic) space- 
time which admits a hypersurface orthogonal timelike Killing vector field, t M , 
whose orbits are complete. Let be a scalar field on (M,g^ v ) which satisfies the 
Klein-Gordon equation 

V%0-m^ = O, (1) 

with V M being the derivative operator associated with the spacetime metric g^. 
Since (M, g^) is static, the equation of motion (1) can be rewritten as 

£*—M. (2) 

where t denotes the Killing parameter and 

A = -VD^VDi) +m 2 V 2 , (3) 

where V = (—t^t^) 1 ^ 2 and D L is the derivative operator compatible with the 
induced metric defined on a spacelike hypersurface £ orthogonal to the orbits of 

We can view A as an operator (with domain C£°(E)) on the Hilbert space, 
L 2 (E, V~ 1 dY*), of square integrable functions on E with volume measure V^dH, 
where <i£ denotes the natural volume element on £ arising from the induced 
metric. Then, provided that m 2 , ^ 0, A is a positive, symmetric operator 
on L 2 (E, V^dY,). Hence, A always admits at least one positive self-adjoint 
extension — namely, the Friedrichs extension, Ap. Let Ap denote any positive 
self-adjoint extension of A. For any initial data (0o?0o) £ L 2 (E, l /_1 rfE) x 
L 2 (E,\/~ 1 rfE), define t by 

4> t = cos(4 /2 t)0 o + A E 1/2 sin(4 /2 t)0 o . (4) 

Then it was shown in [17] that whenever (0 O , 4>o) £ Co°(^) x C^°(E), the quantity 
4>t defines a smooth solution of eq. (1) on all of M. Furthermore, this solution 
agrees with the solution determined by the initial data (c/>o, </>o) within the do- 
main of dependence of initial hypersurface E. In addition, these solutions have 
a conserved, positive energy E = (<p t , Ae^l 2 + (<fit,<fit)L 2 , and satisfy properties 
(iii) and (iv) mentioned above. 

As already stated, the main result proven in [16] was that the above prescrip- 
tion eq. (4) is the only possible way for defining dynamics, subject to require- 
ments (i)-(iv). Hence, we have the following simple algorithm to systematically 
determine the possible boundary conditions at infinity for a scalar field in AdS 
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spacetime: (a) Write the equation of motion in the form (2), thereby identifying 
the operator A, eq. (3). (b) Determine all of the positive, self-adjoint extensions, 
Ae, of A. (c) Interpret the dynamics defined by eq. (4) in terms of boundary 
conditions at infinity. In the present paper, we shall carry out this procedure to 
completion. 

The only property of A used in the above analysis is that it is a positive, sym- 
metric, elliptic operator on a suitable L 2 Hilbert space. Therefore, the analysis of 
scalar field dynamics given in [17] and [16] can be generalised straightforwardly 
to second order, linear, hyperbolic equations satisfied by arbitrary tensor fields 
provided that in static spacetimes 3 these equations can still be put in the general 
form of eq. (2), with A a positive, symmetric, elliptic operator. In particular, it 
can be applied to Maxwell's equations and the linearised Einstein equation pro- 
vided that these equations have been put in the form eq. (2). However, it is not 
straightforward to put Maxwell's equations or the linearised Einstein equation in 
the required form. For example, although Maxwell's equations for the vector po- 
tential, Ap, in the Lorentz gauge have the desired second order, hyperbolic form, 
the "inner product" needed to make the corresponding A symmetric is not posi- 
tive definite, since it yields a negative contribution from the time component, A , 
of Ap. In Minkowski spacetime, one may set A = by a further gauge choice, 
but such a choice is not possible in a general curved spacetime. Similar remarks 
apply to the linearised Einstein equation in the transverse traceless gauge. 

In this paper, we shall analyse Maxwell's equations and the linearised Ein- 
stein equation in AdS spacetime by making use of the spherical symmetry of AdS 
spacetime to decompose the vector potential and metric perturbation into their 
"scalar", "vector", and "tensor" parts with respect to the rotation group. We 
will then expand each part in the appropriate spherical harmonics and will reduce 
the resulting equations to a collection of decoupled scalar wave equations in two 
dimensions, each of which is of the general form (2). In fact, we will show that 
each spherical harmonic component of the scalar wave equation (with arbitrary 
mass and curvature coupling), Maxwell's equation, and the linearised Einstein 
equation all take exactly the same general form, the only difference being the 
values of two parameters that appear in this equation. Therefore, we may de- 
termine the possible boundary conditions at infinity for scalar, electromagnetic, 
and gravitational perturbations of AdS spacetime by studying the self-adjoint 
extensions of a two parameter family of second order differential operators in one 
variable. 

In our analysis, we shall consider (n + 2)-dimensional AdS spacetime (M, g^) 
for all n ^ 1. (Our results also can be straightforwardly extended to n = for 
the case of a scalar field.) We will work in a globally defined coordinate system 

3 The prescription of [17] has recently been generalised to stationary spacetimes [18]. How- 
ever, no uniqueness results similar to those of [16] are presently known for the stationary case. 
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in which the metric takes the form 



I 2 

d s \ n +2) = — (—dt 2 + dx 2 + cos 2 x n fijdz t dzP^ , (5) 

where ^ij(z)dz l dz^ is the metric of the n-dimensional round unit sphere S n and 
£ denotes the curvature radius. The range of x is (0, 7r/2] and the conformal 
boundary is located at x — 0. The coordinate vector field d/dt is a globally 
timelike, hypersurface orthogonal Killing vector field 4 . By use of the spherical 
harmonic decomposition described in the previous paragraph, we will reduce the 
scalar wave equation, Maxwell's equation, and the linearised Einstein equation 
to wave equations in a 2-dimensional AdS spacetime 

ds 2 = (-dt 2 + dx 2 ) , (6) 

sin x 

which is just the space of orbits (0 2 ,g a b) of (n + 2) -dimensional AdS spacetime 
under the rotation group SO{n + 1). 

In the next section, we will perform the spherical harmonic decomposition of 
the equations of motion for scalar fields, electromagnetic fields and gravitational 
perturbations and show that they can be written in the form of a Klein-Gordon 
equation in 2-dimensional AdS spacetime. The analysis of self-adjoint extensions 
of the operator A appearing in this wave equation and the relation between self- 
adjoint extensions and boundary conditions at infinity are given in Section 3. 



Notation and conventions 

We will generally follow the notation and conventions of [6] . However, in this pa- 
per we will make frequent use of the above noted fact that (n+2)-dimensional AdS 
spacetime, (M, g^), can be expressed as a warped product of an n-dimensional 
round unit sphere, (S n ,%j), and a 2-dimensional AdS spacetime, (0 2 ,g a b)- To 
distinguish between tensors on these different spaces, we will use greek indices 
to denote tensors on M, we will use latin indices in the range a,b, ... ,h to de- 
note tensors on O 2 , and we will use latin indices in the range ... ,p to denote 

4 The analysis of dynamics in this paper will be performed with respect to this globally 
timelike Killing field. There are other hypersurface orthogonal Killing fields in AdS spacetime 
that are timelike in a region that includes at least a portion of infinity. We believe that an 
analysis of dynamics with respect to these static Killing fields would give rise to precisely the 
same choices of boundary conditions, but we have not fully investigated this issue. 

There have arisen some interests in implications of the topology for the AdS-CFT correspon- 
dence [19, 20, 21]. The simplest situation can be obtained by considering a quotient of AdS 
spacetime by a certain discrete subgroup of the group of the AdS isometries. We believe that 
an analysis of dynamics in such a geometry parallels the present analysis performed in the 
universal covering space. 
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tensors on S n . When convenient, we will also use these indices to denote coor- 
dinates and coordinate components in these spaces. Thus, we may occasionally 
write equations like eq. (5) in the form 

g^dxV = g ab (y)dy a dy b + r^y^z^dz'dz^ , (7) 

where r = icosx/ sin a;. 

We will denote the metric compatible derivative operators on (M, g^ u ), (S n , 7^), 
and (0 2 ,g a b) by V M , Di, and V a , respectively. In order to avoid introducing ad- 
ditional notation, we will also use 7^ and Di to denote the metric and derivative 
operator of the general Riemannian manifolds considered in the next section. In 
a few instances, we will need to introduce the derivative operator associated with 
the hypersurfaces, S, orthogonal to the static Killing field. In order to avoid 
introducing additional notation, we will also use Di to denote this derivative 
operator (as we already did in eq. (3) above). 

2 Field equations in anti-de Sitter spacetime 

In this section, we will analyse the equations of motion for scalar, electromag- 
netic, and gravitational perturbations of AdS spacetime. Our analysis will be 
based on the decomposition of the modes of these fields into their scalar, vector, 
and tensor parts as well as the expansion of these parts in appropriate spheri- 
cal harmonics. The desired general decomposition theorems are derived in the 
first subsection, and the spherical harmonic expansion is discussed in the second 
subsection. The explicit form of the scalar, electromagnetic, and gravitational 
perturbation equations in AdS spacetime are then derived in the final subsection. 
Most of the results of subsections 2.1 and 2.2 are well known and are given in 
many references (see, e.g., [22, 23, 24]), but complete derivations/proofs have not 
usually been provided in these other references. 

2.1 Decomposition of Vectors and Symmetric Tensors on 
Compact Manifolds 

We begin by recalling some general results on elliptic operators on compact man- 
ifolds. Let (A4,7y) be a compact Riemannian manifold. Let P be a differential 
operator of order m (with smooth coefficients) mapping tensor fields of type (k, I) 
on M. to tensor fields of the same type on M.. Thus, P has the explicit form 

P = F" mi ., mk ni - ni Pl ~ Pm D Pl ...D Pm + lower order terms (8) 

where Di denotes the derivative operator associated with 7^. For each x G M. 
and each covector fcj at x, the symbol, &p(x, k), of P is defined as the following 
map from tensors of type (k, I) at x to tensors of type {k, I) at x 

nn ( T U\ — P^l-H . . n 1 ...n lPl ...p m i 7 /q\ 

upyjj, 1%) — 1 n—3i m\...m k r»p! ■ ■ ■ i^p m ■ 
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The operator P is said to be elliptic if the map a P (x, k) is invertible for all x E M 
and all fcj ^ 0. A well known result (usually referred to as "elliptic regularity") 
states the following: Let P be elliptic and let j be a smooth tensor field of type 
(k,l) (where here and in the following we will generally suppress all indices). Let 
u be a distributional tensor field of type (k, I) which satisfies the equation 

Pu = j (10) 

in the distributional sense. Then u is smooth. 

We can use the metric, 7^ to define a natural L 2 inner product on (complex) 
tensor fields of type (k, I) on M. via 

("•«•)/■ / »"'-V.. J ,«n.V 1 - J '. (11) 

Jm 

where the star denotes complex conjugation. (In this equation indices are raised 
and lowered with 7^ and the integration is with respect to the volume element 
associated with 7^.) We may view any differential operator, P, as an operator 
on L 2 , defined on the dense domain of smooth tensor fields. Now let P be an 
elliptic operator on Ai that is symmetric when viewed as an operator on L 2 , i.e., 
for all smooth u and w 

(u,Pw) = (Pu,w) . (12) 

As will be discussed further in subsection 3.1 below, the self-adjoint extensions 
of P are determined by the "deficiency subspaces", i.e., the vectors z± in L 2 that 
satisfy 

P*z± = ±iz± , (13) 

where P* denotes the adjoint of P. However, eq. (13) implies that z± satisfies 
the equation (P^fi)z± = in the distributional sense, so by elliptic regularity z± 
must be smooth. Hence, z± lies in the domain of P, and we may replace P* by P 
in eq. (13). Taking the inner product of this equation with z±, we immediately 
obtain a contradiction unless z± — 0. Consequently, the deficiency subspaces 
are always trivial, and any symmetric, elliptic operator on a compact manifold is 
always essentially self-adjoint, i.e., it has a unique self-adjoint extension. Thus, 
the usually difficult analysis of determining whether a symmetric operator admits 
self-adjoint extensions and of characterising the self-adjoint extensions that do 
exist (see section 3 below) is always trivial for elliptic operators on a compact 
manifold. 

Now let P be any symmetric, elliptic operator on M. that is bounded from 
below as an operator on L 2 . By the previous remarks, P has a unique self- 
adjoint extension, which we shall also denote by P. The operator exp(— sP) can 
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be then denned by the spectral theorem, and, for s > 0, it is a bounded self- 
adjoint operator. Since Pexp(— sP) also is bounded operator on L 2 , it follows 
that exp(— sP) also is a bounded map from L 2 into the mth Sobolev space, H m 
of A4, where m is the order of P. Since the embedding of H m into L 2 is a 
compact map when M. is compact (see, e.g., [25]), it follows that exp(— sP) is 
a compact, self-adjoint operator on L 2 . Hence, by the Hilbert-Schmidt theorem 
[26], it admits an orthonormal basis of eigenvectors, with eigenvalues converging 
to zero. It follows immediately that P itself has a purely discrete spectrum, and 
that all of the eigensubspaces of P are finite dimensional. Note also that by 
elliptic regularity, all eigenvectors of P are smooth tensor fields. Consequently, 
an equation of the form (10) will admit a solution if and only if j is orthogonal 
to kerP, where kerP is the subspace of solutions to Pz = 0. Furthermore, if a 
solution to (10) exists, it obviously is unique up to addition of vectors in kerP. 
As already stated above, we have dim(kerP) < oo. 

We now are ready to state our two decomposition results. The first one is 
essentially a simple case of the Hodge decomposition theorem for differential 
forms. 

Proposition 2.1: Let (At,7y) be a compact Riemannian manifold. Then any 
dual vector field, Vi on M. can be uniquely decomposed as 

Vi = Vi + DiS (14) 

where D % Vi = 0. We refer to Vi and S, respectively, as the vector and scalar 
components of v j. 

Proof: If the decomposition exists, then S would have to satisfy 

D i D i S = D i v i . (15) 

Since P = —D l Di is a positive elliptic operator, we can solve this equation for S 
if and only if the "source term" on the right side is orthogonal to kerP. However, 
kerP consists of only constant functions and since j M D l Vi = 0, the source is 
indeed orthogonal to kerP. Hence, we obtain a solution S, which is unique up 
to the addition of a constant. We obtain the desired decomposition by defining 
Vi = Vi — DiS. This decomposition is clearly unique, since the only ambiguity in 
S is the addition of a constant, which does not affect the decomposition. □ 

The second proposition concerns the decomposition of symmetric, second rank 
tensor fields. 

Proposition 2.2: Let (Ad^ij) be a compact Riemannian Einstein space, i.e., 
Rij = c~fij for some constant c. Then any second rank symmetric tensor field, ty 
on M. can be uniquely decomposed as 

t ■ ■ = T 

+D {i V j} 
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+ \DiDj - - llJ D m D m j S + -l l] t m m , (16) 

where £>%■ = 0, T\ = 0, and D% = 0. We refer to % ^, and (S,t m m ), 
respectively, as the tensor-, vector- and scalar-type components of t^. 
Proof: We define 

Jij — tij mi (I'' 7 ) 

Jj = D l J tj , (18) 
J = D i J i = D'EP Jij. (19) 

Consider the equation 

D l Di + ~^-rc) W = J. (20) 

n — 1 / n — 1 

This equation can be solved for W if and only if the source J is orthogonal to 
the solutions to 

D*Di + - At c J C/ = . (21) 
n—1 / 

Taking the derivative of this equation and using Rij = c^j, we find 

D*Di + —^—c) DM = . (22) 
n — 1 J 

Contracting with D^U and integrating over A4, we find that any solution, U, to 
eq. (21) must also satisfy 

DiDjU - ^% j D m D m U = . (23) 

i.e., D l U must be a conformal Killing field. It is easily checked that any solution, 
U, to eq. (23) is automatically orthogonal to J, so eq. (20) always can be solved 
for W. If c ^ 0, integration of eq. (20) yields f M W = 0. If c = 0, then 
U = const solves (21), so we can always choose the solution W to eq. (20) to 
satisfy J M W = 0. Therefore, in all cases, we can solve the equation 

D i D i S = W (24) 



to obtain S. 

Next, we consider the equation 



ii — 1 

D % D(iVj) = Jj DjW - cDjS . (25) 

Th 



Again, the left side is a symmetric elliptic operator on dual vector fields, so we 
can solve eq. (25) provided that the right side is orthogonal to the kernel of this 
operator. However, it is easily seen that the kernel consists precisely of Killing 
vector fields and that the right side is orthogonal to any Killing vector field. Thus, 
we can always solve eq. (25) to obtain Vj. Furthermore, taking the divergence of 
eq. (25), we find 

n — 1 

DW'D^ = J DWjW - cDWjS = (26) 

Th 

and hence 

(DWi + c) D j Vj = . (27) 

Taking the gradient of this equation, we find 

{D i D i )D m (D j Vj) = (28) 

which implies that DiD m D^Vj = 0. In particular, we have D l Di{D^Vj) = 0, 
which implies that D^Vj is constant. However, since J M D^Vj = 0, we obtain the 
desired result that D^Vj = 0. 

Having now defined S and Vj, we simply define by eq. (16). It is then easily 
verified that T\ = and that D l Tij = 0. Finally, uniqueness of the decomposition 
follows from the fact that each of the 4 terms in eq. (16) are manifestly orthogonal 
in the L 2 inner product. Consequently, if there were two distinct ways of writing 
the decomposition (16) of tij, their difference would yield a formula for the zero- 
vector, expressing it as the sum of 4 orthogonal vectors, not all of which are zero. 
□ 

Finally, we note that compactness of M. was used in the arguments of this 
subsection in the following two essential ways: (1) to argue that any symmetric 
elliptic operator is essentially self-adjoint, and (2) to argue that the spectrum of 
any self-adjoint elliptic operator is discrete. The symmetric, elliptic operators 
considered in the above propositions will be essentially self-adjoint on many non- 
compact manifolds of interest, but their spectra typically will be continuous. 
Consequently, although analogous decomposition results should hold for many 
non-compact manifolds, it will not be possible to give such a simple, general 
proof, i.e., the decomposition results for non-compact manifolds must be analysed 
on a case-by-case basis. 

2.2 Spherical Harmonic Decomposition on an n-Sphere 

Proposition 2.1 of the previous subsection holds for arbitrary compact Rieman- 
nian manifolds, and Proposition 2.2 holds for arbitrary compact Riemannian 
Einstein spaces. However, on general spaces, an operator of interest may "mix" 
the various components of vectors and tensors, e.g., when acting on the purely 
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"tensor part" of a symmetric rank-two tensor field, the resulting tensor field could 
have scalar and vector components. However, we shall see below that such "mix- 
ing" cannot occur on the (round) n-sphere for any operator that is rotationally 
invariant (and, in the case of the 2-sphere, also parity invariant). Consequently, 
the decompositions of the previous subsection are particularly useful in this case, 
since one can then analyse the scalar, vector, and tensor components separately. 
The purpose of this subsection is to establish this result as well as to introduce 
the spherical harmonic expansion of each component. 

Let be the round metric on the n-sphere, S n . Then (S n , 7^) is an Einstein 
space with c = n — 1, so, in particular, Proposition 2.2 applies. Consider, first, 
the scalar spherical harmonics. The rotation group SO(n + 1) acts naturally 
on the Hilbert space L| of square integrable functions / : S n — > C. Since the 
Laplacian, D l Di, is invariant under rotations, each eigensubspace of D l Di must 
also be invariant under rotations. The scalar spherical harmonics, S ks are defined 
to be an orthonormal basis of eigenvectors of D l D h i.e., 

(L>*A + £;|)§ ks = 0, (29) 

and 

/ § ks § k ' s = <5k S) k' s • (30) 
Js n 

An explicit choice of § ks can De found in e.g., [27] and references therein. The 
eigenvalues k s appearing in eq. (29) are k s = 1(1 + n — 1), I — 0, 1, 2, . . .. 

The rotation group also acts naturally on the Hilbert space of square inte- 
grable vector fields on S n . In addition, it is easily seen that rotations map the 
Hilbert subspace, L v , of divergence- free vector fields into itself. Furthermore, 
it is easily checked that the Laplacian, D l Di, also maps L v into itself. Vector 
spherical harmonics, V kv j, are defined to be an orthonormal basis of eigenvectors 
of D l Di in L v , i.e., they satisfy 

(LPDj + k v )V kvi = , D'V^i = , (31) 

and 

/ V kvt V kV l = <J kV)kV . (32) 

An explicit choice of V kv j can be found in e.g., [27] and references therein. The 

eigenvalues k v appearing in eq. (32) are k\ — 1(1 + n — 1) — 1, I — 1, 2, The 

number of the independent components of V kyi is (n — 1), so vector spherical 
harmonics are nontrivial only when n ^ 2. 

Finally, the rotation group also acts naturally on the Hilbert space of square 
integrable rank-two symmetric tensor fields on S n . In addition, it is easily seen 
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that it maps the Hilbert subspace, L|, of trace-free, divergence-free symmetric 
tensor fields into itself. Furthermore, it is easily checked that the Laplacian, 
D l D{, also maps L\ into itself. Tensor spherical harmonics, Tk T y are defined to 
be an orthonormal basis of eigenvectors of D l Di in L 2 ^, i.e., they satisfy 

(D m D m + k 2 T )T kTt:j = 0, ZVV* = 0, IV, = 0, (33) 

and 

/ T kTij T k , T ij = ,5 kT , k , T . (34) 
Js n 

An explicit choice of Tk T ij can be found in Ref. [27]. The eigenvalues k\ ap- 
pearing in eq. (33) are k\ — 1(1 + n — 1) — 2, I = 2, 3, . . .. The number of the 
independent components of T^ Ti j is (n — 2)(n + l)/2, so, in particular, tensor 
spherical harmonics are nontrivial only when n ^ 3. 

The "total squared angular momentum", J 2 , is the Casimir element of the 
universal enveloping algebra of the Lie algebra of SO(n + l) obtained by summing 
the squares of the elements of an orthonormal basis, {£(«)}, of the Lie algebra. 
For tensor fields on S n , the action of J 2 is obtained by replacing each by the 
Lie derivative with respect to the corresponding Killing vector field on S n . The 
action of J 2 on scalar functions on S n is thereby computed to be simply 

J 2 = -D l Di . (35) 

Similarly, the action of J 2 on dual vector fields and on rank-two symmetric tensor 
fields on S n are, respectively, 

(Jy)i j = ~{D m D m - (n - 1))V , (36) 

and 

(4h M = -(D m D m - 2n)5 i k 5 j l - 2 ll3l kl . (37) 

Consequently, for scalar spherical harmonics of order /, J 2 takes the value /(/ + 
n — 1); for vector spherical harmonics of order /, J 2 takes the value 1(1 + n — 
1) + (n — 2); and for tensor spherical harmonics of order /, J 2 takes the value 
1(1 + n — 1) + 2(n — 1). By inspection, the values of J 2 for scalar, vector, and 
tensor harmonics never are equal, except for the case n = 2, when the values of 
J 2 for scalar and vector spherical harmonics agree for all /. Thus — except in the 
case n = 2 — the scalar, vector, and tensor harmonics always span inequivalent 
representations of the rotation group. This implies that there cannot exist any 
(nonzero) rotationally invariant map between any pair of the spaces L\ and 
Ly and L 2 S — with the exception of the case n = 2, where there could exist a 
rotationally invariant map between L v and L\ (see below). In particular, this 
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means that any rotationally invariant operator on rank-two symmetric tensor 
fields on S n cannot take the tensor component into a nonzero vector or scalar 
component; it cannot take the vector component into a nonzero tensor component 
or (if n > 2) a nonzero scalar component; nor can it take the scalar component 
into a nonzero tensor component or (if n > 2) a non-zero vector component. 

When n = 2 the scalar and vector harmonics span equivalent representations 
of SO(3) and, indeed, the vector harmonics can be chosen to be 

Y kl = eij&Sk , (38) 

where is the (positively oriented) volume element on S 2 associated with 7^. 
Thus, one sees explicitly that there exist nonvanishing S*0 (3)-invariant maps be- 
tween Ly and L 2 S . However, as can be seen from eq. (38), the vector and scalar 
harmonics behave differently under parity, and thus span inequivalent represen- 
tations of 0(3). (The scalar representations of 0(3) are said to be of "polar" (or 
"even parity") type, and the vector representations are said to be of "axial" (or 
"odd parity") type.) Consequently, even when n = 2, no rotationally and parity 
invariant operator can "mix" the scalar and vector components of a vector field 
or symmetric, rank-two tensor field. 



2.3 Scalar fields on AdS Spacetime 

We now turn our attention to the derivation of the desired form of the field equa- 
tions for scalar, electromagnetic, and gravitational perturbations of AdS space- 
time. As previously noted, the metric (5) of (n + 2) -dimensional AdS spacetime 
takes the form of a warped product of an n-dimensional base space (S n ,~fij) and 
a 2-dimensional space of orbits (0 2 ,g a b): 

ds 2 {n+ 2) = g^dx"dx u = g ab (y)dy a dy b + r%) 7 (z)^G^ , (39) 

where 7^ is the round, unit sphere metric on S n , and 

r = I . (40) 

sin a; 

In the coordinates (t, r) the metric of the space of orbits takes the form 

g ab (y)dy a dy b = -V 2 l 2 dt 2 + — , V 2 = 1 + - . (41) 

Our general strategy will be to project all tensor fields into parts that are tangent 
and orthogonal to S n , to decompose all tensor fields on S n into their scalar, vector, 
and tensor components, and to expand each component on S n in the appropriate 
spherical harmonics. 

To avoid confusion when comparing with other analyses, it should be empha- 
sised that our analysis involves an "n + 2" decomposition of AdS spacetime rather 
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than an "(n + 1) + 1" decomposition. Thus, for example, our notion of the scalar, 
vector, and tensor components of a metric perturbation are defined with respect 
to the action of the isometry subgroup SO(n +1), whose orbits are n-spheres. 
It would also be possible to write AdS spacetime as a Robertson- Walker space- 
time and decompose tensors with respect to the isometry subgroup acting on the 
(n + l)-dimensional surfaces of homogeneity and isotropy in the manner done, 
e.g., in [22, 23]. This would give rise to completely different notions of the scalar, 
vector, and tensor components of a metric perturbation. 

The simplest case is that of a Klein-Gordon scalar field, 0, where no projec- 
tions are necessary. We expand in terms of scalar spherical harmonics §k(-2 l ) 

0(t,x,^)=r-/ 2 ^$ k (t,r)§ k (^). (42) 

The equation of motion for each $k(£, t) is obtained by straightforward substitu- 
tion into (1). Returning to the original coordinates (t,x), we obtain 

** = (JL-?=W-£=1I£\« ,43) 

dt 2 \dx 2 sin 2 x cos 2 x J 

where the index k labelling the harmonics has been omitted. Here the constants 
v and a in the effective potential term are given by 

2 1 n{n + 2) 22 
v 2 - - = - + m\i 2 , (44) 

and 



2 1 in , i\ , n(n-2) 



<j*-- = l(! + n-l) + -^- r ^. (45) 

Since the scalar curvature, R = — (n+ l)(n + 2)£~ 2 , is constant in AdS space- 
time, the modification to eq. (1) resulting from the addition of the curvature 
coupling term — £R is equivalent to changing the mass by itlq — > m 2 , where the 
effective mass, m 2 , is defined by 

m 2 = ml - £(n + 1) (n + 2)£~ 2 . (46) 

Note that for the case of a conformally coupled scalar field (i.e., m\ = and 
£ = n/4(n + 1)), we have 

J = \- (47) 



14 



2.4 Electromagnetic fields 

Maxwell's equations for the electromagnetic field tensor F^ v are 

V v F Vfi = , (48) 
V [CT F H =0. (49) 

Since AdS spacetime is topologically trivial the second equation (49) immediately 
implies the existence of a vector potential A^, so that 

= V»A V - V V A„ . (50) 

Equation (50) is the general solution to eq. (49). 

We may write A^ as the sum of its projection orthogonal to S n and its pro- 
jection tangent to S n . The components of A^ orthogonal to S n behave as scalars 
under rotations of S n . The projection of A^ tangent to S n can then be de- 
composed into its scalar and vector parts in accordance with Proposition 2.1. 
Maxwell's equations are invariant under rotations and parity transformations, 
and hence cannot couple the vector and scalar parts of A^. Therefore, we may 
treat these cases separately. 

In order to express Maxwell equation (48) explicitly, it is useful to write down 
the non- vanishing components of the Christoffel symbol Y^ v \ associated with g^ u , 

r a fe c = r a f> c , T a ij = —rV a r^ij , T\j = ——S l j , F l jk = P jfe , (51) 

where T a i JC and f ^ are the components of the Christoffel symbols associated 
with the metrics g a b and 7^, respectively. 

Vector-type component 

The purely vector part, A^, of A^ can be expanded in vector spherical harmonics 

as 

Aldx» = J2<Pv^dz\ (52) 

Since the gauge freedom in A^ is purely scalar, it follows that A^ — and hence 
each 0y k — is gauge- invariant. Substituting eq. (52) into eq. (48), we obtain 

-Lv c (r»-^,M - ^ + <rV = , (53) 

where we have omitted the index k here and in the following. In terms of the 
rescaled variable $y = A n ~ 2 ^ 2 (j)v, we obtain 

V»V A- - 1 + Izl^lM 4,, _ 4 + = . (54) 

2 [ r 2 r z J r z 
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We can easily check that the a-components of eq. (48) provide just a trivial 
equation, as they must since they are scalar components of Maxwell's equations. 

Substituting V a V a r/r = 2/f and (V c r)V c r/r 2 = l/i 2 + 1/r 2 into eq. (54) 
we obtain 



n(n — 2) 



1 



$ y = 0. (55) 



In terms of the original coordinates (t,x), eq. (55) takes exactly the same form 
as eq. (43), with 



n(n — 2) 1 



^ = ^^ + 4' (56) 
and with a 2 again given by eq. (45). 

Scalar-type component 

The scalar part, A^, of A^, consists of two quantities: a vector, A a , orthogonal 
to S n and the scalar component, A, of a vector tangent to S n (see Proposition 
2.1). We may expand both of these quantities in scalar spherical harmonics §k, 

Aldx» = {A^dy a + AkASktfe*) , (57) 

The contribution of each spherical harmonic component to the field strength is 
given by 

F ab = (V a A b -V b A a )§, (58) 
F ai = (y a A + k s A a )D i S, (59) 
^ = 0, (60) 

where we have again omitted writing the index k. The i-components of eq. (48) 
then take the form 

V a |r n - 2 (v a A + k s A a ^j }=0. (61) 

This immediately implies that there exists a scalar field 4>s such that 

V a <Ps = e ab r n - 2 (V b A + k s A b ), (62) 

where e ab is the metric compatible volume element on (O 2 , g ab ). Since (ps is related 
to the field strength via eq. (59) as 



F a , = — e ab VVsA§, (63) 
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gauge freedom of 4>s is restricted to only the replacement 

05 — ¥ (fis + const . 
In terms of <f>s, the a-components of eq. (48) are written as 



0. 



(64) 



(65) 



So, we have 



r"" 2 V a 



„n— 2 



(66) 



Since the integration constant c in the right hand side of this equation can al- 
ways be absorbed in <f> s by using the remaining gauge freedom (64), we have the 
following equation of motion for the now gauge-invariant variable <p s , 



r n - 2 V a 



In terms of the rescaled variable $5 = r~( n ~ 2 ^ 2 <j)s, we have 



kj 

.,2 



V> 5 ) -^05 = 0. 



(67) 



V a V a $ 5 + 



n - 2 J V a V a r n (V c r)V c r 



fc2 

$ 5 - -f $ 5 



(68) 



Again substituting V a V„r/r = 2/£ 2 and (V c r)V c r/r 2 = l/£ 2 + 1/r 2 , we obtain 



„ T , (n-2)(n-4) 1 



n(n — 2) 



+ Z(Z + ra- 1) 



$5 = 0. (69) 



In terms of the original coordinate system (t,x), this equation again takes the 
canonical form (43) with 



(n-2)(n-4) 1 
4 4 ' 



(70) 



and a 2 again given by eq. (45). Notice that in 4-dimensional AdS spacetime (i.e., 
n — 2), both vector- and scalar- type components of the electromagnetic field 
satisfy the same equation as a conformally invariant scalar field. 



2.5 Gravitational perturbations 

We shall basically follow the derivation of equations of motion for the master 
variables of gravitational perturbations developed in Ref. [28]. The resultant 
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equations are equivalent to those given in [29] in maximally symmetric back- 
grounds and those in [24, 30] in more general backgrounds. 

Gravitational perturbations are described by a symmetric tensor field = 
5g^ u on the background AdS spacetime. We may project relative to S n as 



doc doc 



h ab dy a dy b + 2h ai dy a dz l + hijdz l dz j . 



(71) 



The projection, h a i>, orthogonal to S n in both indices is purely scalar with respect 
to rotations. The orthogonal-tangent projection, h ai , and the tangent-tangent 
projection can be further decomposed into their scalar, vector, and tensor 
parts as 

h ai = D t h a + h^ ai , (72) 



h 



n 



with 



+h Lllj + \D i D j - \ijD m D r , 



£>W 2 >„ = = , 



h 



(o) 



= . 



(73) 

(74) 

(75) 

(i) 



Thus, the tensor part of is h^ 2 \j, the vector part of h^ v consists of (h^ ai , h T j), 
and the scalar component of h^ v consists of (h ab ,h,L,h^)- Each tensorial com- 
ponent decouples from the others and can be expanded in the appropriate tensor 
harmonics, namely Tk T «j, Vk v i, and §>k s > respectively. 



Tensor-type perturbation 

We expand the tensor-type perturbation in spherical harmonics as, 



5> 



(2) 
Tk 



(76) 



k T 



with H^p being a function on (0 2 ,g a b)- Since the gauge freedom 

h^u -> - - , (77) 

does not contain any tensor part, it is clear that H^p is gauge-invariant. Writing 
we find that the linearised Einstein equations yield 



V a V a $ T - 



n(n + 2) 1 



n(n 



+ l(l + n-l) 



4 I 2 

This equation again takes the form (43) with v 2 given by 

1 



$ T = 0. 



2 = n(n + 2) 

4 4' 



(78) 



(79) 



and a 2 again given by eq. (45). 
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Vector-type perturbation 



Vector-type components can be expanded in terms of vector spherical harmonics 

as 



1,(1) . _ .v t - 

ll - ai — / - t - t k a v kj ' 



(1) 
T j 



ki • 



(80) 
(81) 



Note that, as shown in the Appendix, I = 1 modes do not correspond to 
dynamical perturbations, so hereafter we assume I ^ 2. 

The vector component, ^Yidz\ of the gauge transformation eq. (77) yields 



if, 



(i) 



4 1 ) - 



(82) 
(83) 



It follows immediately that a natural gauge-invariant combination is given by 



Z n = ff « - r 2 V, 



ff 



(i)' 



(84) 



Substituting this into the linearised Einstein equations, we have 

V„ (r"- 2 Z") = , (85) 

k *'- {n - l) -Z„ = . (86) 



„n+2 



These equations have the same form as the 2-components and a-components, 
respectively, of eq. (48) in the case of the scalar part of electromagnetic fields, 
with V a A + k s A a replaced by Z a . In parallel with the introduction of (ps in the 
electromagnetic case (see eq. (61)), we find from eq. (85), that there exists a 
scalar potential (py such that 



Z a = -^ e abV b (pv ■ 



Substituting this into equation (86), we obtain 

1 



r n+1 V b [-V b 4>v)-{k v -{n-l)}<t>v 



0. 



(87) 



(88) 



From this equation (after eliminating an integration constant using the degree of 
freedom 0y — > (fiy + const, as we have done in the electromagnetic field case), we 
obtain 



V a V a $y - 



n(n-2) 1 

4 £2 



+ 



n(n - 2) 



+ l(l + n-l) 



1 



$y = , (89) 
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where $y = r n / 2 <f>v- This equation again has the form (43) with v 2 given by 



n(n — 2) 1 



^ = "V^ + 4' (90) 

and a 2 again given by eq. (45). Note that this equation is precisely the same as 
the vector- type of Maxwell's equations (55). 

Scalar-type perturbation 

The scalar components of h^ v can be expanded in scalar spherical harmonics as 
follows: 

hat = ^ifJVSk, (91) 
h a = ^^ 0) a -S k , (92) 

h L = E^S-Sk, (93) 

k S 

4 0) = E^'Sk, (94) 

k S 

where each of the expansion coefficients, (# (0 U, # (0) a , i^| 0) , ^ 0) ) are functions 
on (0 2 ,g ab ). 

Note that the I = = 0) mode preserves the spherical symmetry of the 
background and, by Birkhoff 's theorem, cannot describe any dynamical degrees of 
freedom. For I = 1 (k§ = n), the corresponding perturbations also do not describe 
any physical dynamical degrees of freedom, as we show in the Appendix. Hence, 
hereafter we consider only the / ^ 2 modes. 

The scalar-type gauge transformation is generated by ^ a Eidy a + ^D i 'Sdz\ 
Observing that the scalar-type components transform as 

H^ ab H^ ab -V a ^\-V b ^ a , (95) 
H (0 \ - H^ a - e\ ~ r 2 V a (^pj , (96) 

H f> -> i/( ) + ^(o)-2r(VV)e (0) a , (97) 
- H®-2g°\ (98) 
one can check that the following combinations are gauge-invariant: 

Z = nr n ~ 4 j#i 0) + + 2r(V a r)X a \ , (99) 



" n " 2 (0 U + V a X b + V 5 X a ) + y -^-A Z g ah , (100) 
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where 



X a ee -H^ a + ir 2 V fl ) , (101) 

which itself is gauge-dependent: X a — > X a + £^ a . 

The full set of linearised Einstein equations for the scalar-type perturba- 
tions are written out in eqs. (63) — (66) of Ref. [28], where one must substitute 
F = r~ n+2 Z/(2n) and F ab = r' n+2 Z ah - {(n - 1) jn}r~ n+2 Zg ah in order to ex- 
press these equations in terms of the variables we have introduced here. The 
(a, i)-components and the traceless part of the (i, j)-components of the linearised 
Einstein equations yield 

V b Z\ = V a Z c c , (102) 
Z\ = Z. (103) 

We will now show that the variable Z a b (and, hence, by eq. (103), also the 
variable Z) can be replaced by a single scalar field, 0£. To do so, we will want 
05 to satisfy 

V"Va " |) 05 = Z\ (104) 

on the space of orbits (0 2 ,g ab ). However, (0 2 ,g ab ), of course, fails to be globally 
hyperbolic, and although the methods of this paper allow us to prove existence 
of solutions to the homogeneous equation with initial data in L 2 (S, V^gE), we 
cannot guarantee existence of solutions to the inhomogeneous equation with an 
arbitrary smooth source that might, in particular, fail to lie in L 2 (£,K -1 <i£) 
at each time. Nevertheless, we can overcome this difficulty by noticing that 
(O 2 , — gab) is a globally hyperbolic spacetime, as can be seen by inspection of 
eq. (6), recalling that the coordinate ranges of (x,t) are — oo < t < oo and 
< x ^ 7r/2. Thus, we can solve eq. (104) for arbitrary smooth Z a a by specifying 
arbitrary smooth "initial data" (4>so(t), (d(f>so/dx)(t)) on a hypersurface x = xq. 
Furthermore, we will specify this initial data so as to satisfy 

t h (v a Vb - ^g a ?j 0s = t b Z ab (105) 

on the initial data surface x = xq, where t a denotes the static Killing field (d/dt) a . 
(The above equation yields a system of two linear ODE's in t for <pso arid d<pso/dx 
and hence always can be solved.) Now define 

Tab = Z ah ~ (VaV b ~ ^g a b) 05 ■ (106) 
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Then T ab is transverse and traceless on (0 2 ,g ab ), i.e., 

V 6 T 6 a = 0, T\ = 0. (107) 

We define 

v a = T ab t b . (108) 

It follows immediately from the Killing's equation V ' a h + Vf>t a = and eq. (107) 
that 

V a v a = 0, (109) 

and, consequently, there exists a scalar field s such that 

v a = e ab V b s. (110) 

By our above choice of initial data for <p s , we have v a = at x = x , so V b s = 
at X — Xq . By adding a constant to s if necessary, we therefore can ensure that 
on the "initial data surface" x = x we have s = and ds/dx = 0. 
On the other hand, we also have 

V a (e ab v b ) = V a (e ab T b %) 
= V a (e cb T b a t c ) 
= e cb T b a V a t c 

= —-^ ch T b a e ac e de V d t e 
= U b a T b a e de V d t e 

= 0, (111) 
where we have used the fact that for any traceless tensor T ab on (0 2 ,g ab ), 

T a b e bc = T c b e ba (112) 

at the second line and the formula V a t c = — (1 /2)e ac e de \/ dt e for any Killing vector 
field t a on (0 2 ,g ab ) at the fourth line. Substituting eq. (110) into eq. (Ill), we 
obtain 

V a V a s = 0. (113) 

Consequently, by Cauchy evolution in x, we obtain s = and, hence, v a — 
throughout (0 2 ,g ab ). However, using the formulae 

g ab = ^(-tah + t ac t bd t c t d ) , (114) 
£ a b£cdt d = t a g bc — t b g ca , (115) 
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we find 



Tab — T a c g, 



cb 



— (—T a c t c tb + T d c e ca eb e t d t e ) 

= ( Vatb + Vbta ~ 9abVct C ) 

= 0, (116) 

where we have used eq. (114) at the second line and eq. (115) at the third line. 
This proves that we can express the gauge-invariant perturbation variable Z ab as 



J ab 



V a V 6 - 4?a*>Us- (117) 



Note that the residual freedom of 0s, which leaves Z ao unchanged, is 

05^05 + 00, (v a V 6 - ^g ab ^j 0o = . (118) 

Substituting the expression (117) and eq. (103) into the (a, 6)-components of 
the linearised Einstein equations given in Ref. [28], we obtain 

(v a V 6 -^ a ^£(0s) = O, (119) 



where 



E{<p s ) ee r 2 ( V a V a - n^V a - + ^ | 5 . (120) 



r 



Viewing eq. (119) as three differential equations for E, we find the general solu- 
tion, 

E S oi = Cir + V (c 2 sin t + c 3 cos t) , (121) 



where ci, C2, C3 are arbitrary constants and V = (—tHa) 1 ^ 2 = \Jl + r 2 /£ 2 . 

On the other hand, since eq. (118) is the same as eq. (119), the residual 
freedom O is also given by 

O = c [r + V (c' 2 sint + c' 3 cost) , (122) 

with arbitrary constants d x ,d 2 ,d z . Using these explicit forms of E so i and O , we 
see that under the change (118), E(4> s ) transforms as 

E(<j> S ) -> ^(05 + 00) 

T 2 ( T 2 \ 

= E so i - (k 2 s - n)0 o + n—(f) - nr fl + — J <9 r o 
= £ soi - (fc§ - n)0 o - nqr 

= c"r + V (4 sin * + eg cos i) , (123) 
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where c'[ — c± — kgC^ = c 2 — (k s — n)d 2 c'l = c 3 — (k 2 s — n)c' 3 . Thus, setting these 
arbitrary constants c", 0%, C3 zero, we can make E(<f)s + 0o) vanish. Therefore, 
from eq. (120), we obtain 



V a V a $ s 



(n-2)(n-4) 1 



+ 



mn 



2) 



+ l(l + n-l) 



$ 5 = 0, (124) 



where $5 = r~ n/2 (0s + O )- 

Again, the equation for scalar-type gravitational perturbations takes the form (43) 
with v 2 given by 



= ( w -2)(n-4) 1 
4 4 ' 



(125) 



and a 2 again given by eq. (45). 

Thus, we have reduced the problem of solving the gravitational perturbation 
equations to solving the 3 scalar equations (78), (89), and (124) for <E>r, and 
$5, respectively. Note that the vector- and scalar-type gravitational perturba- 
tions in 4-dimensions (i.e., n = 2), and scalar-type perturbations in 6-dimensions 
(i.e., n — 4) satisfy the same equation as a conformally invariant scalar field, 
while in all dimensions the tensor-type gravitational perturbations satisfy the 
same equation as a minimally coupled massless scalar field. 

Finally, we note that the tensor spherical harmonics have (n — 2)(n + l)/2 
independent components, the vector spherical harmonics have (n— 1) independent 
components, and the scalar spherical harmonics have 1 independent component. 
Thus, the total number of independent components is (n + 2)(n — l)/2, which 
corresponds to the number of dynamical degrees of freedom for gravitational 
radiation in an (n + 2)-dimensional spacetime. 



3 Dynamics in anti-de Sitter spacetime 

We have seen that the field equations for scalar fields, electromagnetic fields and 
gravitational perturbations can be reduced to a wave equation of the form 

d 2 ( d 2 v 2 -1/4 a 2 -1/4 
dt 2 \dx 2 sin 2 x cos 2 x 

which is defined on the interval (0, vr/2), where x = tt/2 corresponds to the origin 
of spherical coordinates and x = corresponds to infinity in AdS spacetime. In 
all cases, a takes the value (see eq. (45)) 

71 — 1 

(7 = 1 + — . (127) 

On the other hand, the value of v 2 depends upon the field and the type (i.e., scalar, 
vector, or tensor) of perturbation. For scalar fields, u 2 can take any (positive or 
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negative) value, depending on the value of the effective mass, m 2 = m„ + £R, of 
the scalar field. For the electromagnetic and gravitational cases, the values of v 2 
are given by eqs. (56), (70), (79), (90), and (125) above. In this section, we will 
determine all of the possible boundary conditions that can be imposed at infinity 
[x = 0) so as to yield sensible dynamics for this equation. 



3.1 The operator A and its properties 

The general prescription of [17] for defining dynamics (see section 1 above) cor- 
responds to viewing the differential operator 



dx z sin x cos z x 



appearing in eq. (126) as an operator on the Hilbert space Ti = L 2 ([0, tt/2], dx). 
Indeed, taking account of the field redefinition (42), it is easily seen that for each 
spherical harmonic, the inner product for <j)(t, x, z l ) in L 2 (S, V~ 1 dH) corresponds 
to the inner product for $(t,x) defined on H, i.e., 

/ i)*<\> V- x dH = [ 1 dx. (129) 
Jt, Jo 

It is easily seen that, when defined on the natural domain C£°(0, 7r/2), A is a 
symmetric operator. For later purposes, it will be useful to classify A according 
to the possible values of v 2 as follows: 

Case (i): v 2 ^ 1. This case occurs for scalar fields whose effective mass squared 
satisfies m 2 ^ — (n + 3)(n — l)/4£ 2 . In particular, this case encompasses the min- 
imally coupled massless scalar field in all dimensions. This case also encompasses 
the tensor-type of gravitational perturbations in all dimensions where they are 
defined (namely spacetime dimensions greater than or equal to 5, i.e., n ^ 3), 
the vector-type of electromagnetic and gravitational perturbations in spacetime 
dimensions greater than or equal to 5 (i.e., n ^ 3), and the scalar- type of electro- 
magnetic and gravitational perturbations in spacetime dimensions greater than 
or equal to 7 (i.e., n ^ 5). 

Case (ii): < v 2 < 1. This case occurs for scalar fields whose effective mass is in 
the range —(n + 1) 2 /A£ 2 < m 2 < —{n + 3)(n - l)/4£ 2 . (Note that in 4 spacetime 
dimensions, this relation becomes — 9/4£ 2 < m 2 < — 5/4£ 2 , which corresponds to 
the case analysed by Breitenlohner and Freedman [8, 9].) In particular, this case 
encompasses the conformally invariant scalar field in all spacetime dimensions. 
This case also encompasses the vector-type of electromagnetic and gravitational 
perturbations in 4-spacetime dimensions (i.e., n — 2), and the scalar-type of 
electromagnetic and gravitational perturbations in 4- and 6-spacetime dimensions 
(i.e.,n = 2, 4). 
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Case (iii): v 2 = 0. This case occurs for a scalar field whose effective mass 
squared takes the critical value m 2 c = — (n + l) 2 /4£ 2 . It also occurs for scalar- 
type of electromagnetic and gravitational perturbations in 5-dimensional (i.e., 
n — 3) AdS spacetime. 

Case (iv): v 2 < 0. This case occurs only for a scalar field whose effective mass 
satisfies m 2 < —(n + l) 2 /4£ 2 . 

The positivity properties of A are governed by the following proposition: 

Proposition 3.1: If v 2 ^ (i.e., in cases (i), (ii), and (iii) above), then A is a 
positive operator. On the other hand, if v 2 < (i.e., in case (iv) above) then A 
is unbounded below. 

Proof: Consider, first, the case v 2 ^ 0. Let v denote the positive square root of 
v 2 and write 

G„(x) = (cosx) CT+1/2 • (sinx)-" +1/2 . (130) 
We define the operators l) ± by 

D±EE±A_ Gir l^, (131) 

dx dx 

Then A can be expressed as 

A = E>-D + + (-is + cr + l) 2 . (132) 
For arbitrary $ e Dom(A) = C^°(0, vr/2), it is then easily seen that 

($,A$) L2 = (-v + a+ 1) 2 / \$\ 2 dx+ \D + $\ 2 dx, (133) 

Jo Jo 

which establishes the positivity of A. 

Next, we consider the case v 2 < 0. Let / be the C 1 function 

fix) = J( sin;r ) 2 0<x^r ? 

1/1 ; \ (cosrr) 2 {p(smx) 1/2 + q(smx) 1/2 } i]<x^n/2 ' v ; 

where < rj < ir/2 and 

(sin 7]) 3/2 {4 + (cost]) 2 } (sin 77) 5/2 { -4 + (cost/) 2 } 

P= 2{^W ' q= 2(co^ • (135) 

Straightforward calculation shows that for small 77, 

/ ~7TT + • 2 + 2^)f dx / f dx 

\ ox 2 sin x cos^x / /Jo 

= -5z/ 2 logr/ + C + 0(r7), (136) 
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where C is a constant independent of r\. Hence, for v 2 < 0, by taking 77 — > 0, 
the right hand side of eq. (136) can be made to take arbitrarily large nega- 
tive values. However, this does not yet prove the desired result, since / ^ 
Dom(A) = C£°(0, 7r/2). To remedy this, we smooth / as follows. Let x be 
an arbitrary smooth function with support contained within (— 7r/4, 7r/4) satisfy- 
ing J x(x)dx = 1. Define Xe(^) — e" 1 ^^/ 6 ) ( so that Xe approaches a 5-function 
as e — > 0) and define 

fe{x) = f Xe(x - x')f(-e + (1 + e)x')dx' , (137) 

where we have extended the domain of definition of / beyond [0, 7r/2] by setting 
f(x) = whenever x [0,7r/2]. Then, for sufficiently small e, we have f e G 
C^°(0,7r/2) = Dom(A), so formula (133) holds. However, since / is C 1 and goes 
to zero sufficiently rapidly at x — 0, 7r/2, it is straightforward to check that as 
e — > 0, / e and l) + / e converge uniformly on [0, 7r/2] to / and D + f, respectively. 
It then follows that (/ e , Af e ) L 2/\\f e \\^ 2 converges to the left side of eq. (136) as 
e -> 0. Consequently, there exist * G Dom(A) such that A*) L2 /||^||^ 2 takes 
arbitrarily large negative values. □ 

Thus, in case (iv), no positive, self-adjoint extensions of A exist, and there 
is no way of defining a sensible, stable dynamics. On the other hand, at least 
one positive, self-adjoint extension (namely, the Friedrichs extension) must exist 
in cases (i), (ii), and (iii). Our task is to find all of the positive, self-adjoint 
extensions of A in these cases. 

There is a well known systematic method, due to von Neumann, for obtaining 
all of the self-adjoint extensions of an arbitrary symmetric operator on a Hilbert 
space 7i. We will now review this method (see e.g., Ref. [26] for proofs and 
further details). Let A be a symmetric operator with domain Dom(A), and let 
A* denote the adjoint of A. Consider the subspaces J^± C H (called deficiency 
subspaces) spanned by $± G Dom(A*) which satisfy respectively 

A*$ ± = ±i$ ± . (138) 

Let n± = dim(j£±). (Either or both of n + and n_ may be infinite.) If n + 7^ n_, 
then no self-adjoint extension of A exists. If n + = n_ = 0, then there exists 
a unique self-adjoint extension, obtained by taking the closure of A (see e.g., 
Ref. [26]). On the other hand, if n + = n_ 7^ 0, then all of the self-adjoint 
extensions of A can be constructed as follows: Let U : Jtf + — > Jff- be a partial 
isometry, i.e., ||$+|| = ||C/$ + || for all $ + G J^ + . Given U, we define the associated 
self-adjoint extension, Ae, of A by first taking the closure, A, of A and then 
extending the domain of A by adding to it all elements of the form &u = $ + + 
i.e., we define 

Dom{A E ) = {$ + $+ + C/$+|$ e Dom(A), $+ G Jf + } . (139) 
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For all $ G Dom(A E ), we define 



A E $ = A% + i$ + - iU§ 



(140) 



Then A E is self-adjoint, and all self-adjoint extensions of A can be obtained in this 
manner. Thus, the self-adjoint extensions of A are in one-to-one correspondence 
with the partial isometries U : J£+ — > J£L. 

From the definition of A*, it can be seen immediately that eq. (138) is equiv- 
alent to the statement that for all \I> € Dom(A), we have 



In our case, A is an elliptic differential operator with domain Cq°(0, 7r/2) on the 
Hilbert space 7i = L 2 ([0, vr/2], dx). It follows immediately that $± satisfies the 
equation (A =p i)$-i- = in the distributional sense. By elliptic regularity (see 
subsection 2.1 above), it follows that $± must be smooth on (0,7r/2). Thus, in 
our case, the solutions to eq. (138) are precisely the smooth functions on (0, 7r/2) 
that satisfy the ordinary differential equation 



and are square integrable on (0, 7r/2). We turn now to an analysis of solutions to 
this equation. 

3.2 General solutions to the eigenvalue equation 

In this subsection, we analyse solutions to the ordinary differential equation 



with a given by eq. (127), with v an arbitrary non-negative real number, v ^ 0. 
This encompasses cases (i), (ii), and (iii) of the previous subsection. As discussed 
in the previous subsection, the self-adjoint extensions of A can be determined by 
finding all the square integrable solutions to this equation for uj 2 = ±i. More 
generally since the problem of finding the self-adjoint extensions of A is obviously 
equivalent to the problem of finding the self- adjoint extensions of a A for any 
positive real number a, the self-adjoint extensions of A can be determined by 
finding the square integrable solutions to eq. (143) for uj 2 = ±ai for any a > 0. 
It will turn out to be slightly more convenient to study this equation for a = 2, 
and we shall do so below. However, we also will be interested later in solutions 
to eq. (143) for u 2 < 0, i.e., uj imaginary. Therefore, at this stage we will allow 
uj to be an arbitrary complex number. 



((A±i)tt,$±)„ = 0. 



(141) 




(142) 




(143) 
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The general solution to the equation (143) is given, in terms of hypergeo metric 
functions, by 

$ = (sinx) u+1/2 • (cosx)- a+1/2 

x • (cosx) 2a ■ F (q ff ,C£, 1 + <r; cos 2 x) + B 2 ■ F(cos 2 x) } (144) 

where Bi, B 2 G C, and here and hereafter we use the abbreviation 

If a is not an integer (i.e., when the number of spacetime dimensions is even), 
the second solution, F(cos 2 x), is given by 

F(cos 2 x) = F (C%, C"- CT , 1 - o- cos 2 x) . (146) 

If a + 1 is a natural number (i.e., when the number of spacetime dimensions is 
odd) then — with the exception of the case (which we shall not consider further 
here) where u is real and either < ^ a or < < a — then we have [31] 

F = F (C^, C£> 1 + °\ cos 2 x) • log(cos 2 x) 

+ E (T+^r fc! • {W - ^ (0)} • icosx)2k 

k=l ^ 'k 

-E f i*~w£\ - (cos " r2i - (147) 



where 



(0. - (148) 

i/>(z) = ^ogT(z), (149) 
h(k) = il>(£ t<T + k)+1>({^ + k)-1>(l + a + k)-1>(k + l), (150) 

and for the a = case, the third term in the right side of eq. (147) should be 
ignored. 

We consider, now, the behaviour of these solutions near x = ir/2, which 
corresponds to the origin of spherical coordinates in AdS spacetime. For a ^ 1, 
it follows immediately from the general form (144) of the solutions that the square 
integrability of $ at x = n/2 requires B 2 = 0. However, the situation is more 
subtle for the cases a = 1/2 and a = 0. 

For the case a = 1/2 — which occurs only for an S-wave (I = 0) for a scalar field 
in 4-dimensional (i.e., n — 2) spacetime — the leading behaviour near x = n/2 is 
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<&(#) = const. Hence $ is square-integrable even when £? 2 7^ 0, and it might 

appear that both solutions to eq. (143) are acceptable. However, this is actually 
an artifact of our having artificially removed the origin from the spacetime when 
we used spherical coordinates to separate variables. (A similar phenomenon oc- 
curs when analysing the spherically symmetric eigenstates of the Hydrogen atom 
in ordinary quantum mechanics.) If we put the origin back and analyse <3> on 
the 3-dimensional space E, we find that the solution is not acceptable when 
B 2 7^ 0. To see this explicitly, define a Cartesian coordinate system, {X,Y,Z}, 
in a neighbourhood of the origin x = tt/2 on a three-dimensional timeslice, E, 
such that tt/2 — x = \JX 2 + Y 2 + Z 2 . Then, from eq. (42), the corresponding 
field, 4>(X, Y, Z), defined on E is expressed in this neighbourhood as 

\ nil D 
sin x \ ' „ .Do 



</>=[ $Sps - (151) 

\cosxJ ^X 2 + Y 2 + Z 2 

This is square integrable in a neighbourhood of the origin X = Y = Z = 
with the integration measure V~ l dY> = dXdYdZ. However, let A denote the 
3-dimensional operator on L 2 (E, V~ 1 dE) corresponding to A. the leading order 
behaviour of A<p near the origin is given by 

-VD\VD4) ~ -(d 2 /0X 2 + d 2 /dY 2 + d 2 /dZ 2 )(j) 

= 4vr • B 2 ■ 5(X)5(Y)5(Z) , (152) 

where S(-) is the Dirac delta function. As a consequence of this additional 5- 
function source, the "solution", 0, on E when B 2 7^ 0, does not actually cor- 
respond to a distributional solution of A(f) = u 2 <f) but rather to an equation 
that contains an additional 5-function source at the origin. In particular, when 
id 2 = ±2i, we do not obtain a solution to the analog of eq. (138) on E. Therefore, 
only the solution with B 2 = is acceptable. 

Similarly, for the case a = — which occurs only for the S-wave (I = 0) for 
a scalar field in the 3-dimensional (i.e., n — 1) case — $ behaves near x = tt/2 
as <3> ~ (cosrr) 1 / 2 {B 2 log(cos 2 x) + B\ + • ■ •}, which is square-integrable near x = 
n/2 even when B 2 7^ 0. However the original field defined on a 2-dimensional 
surface E behaves near the origin as 

( p^B 2 \og(X 2 + Y 2 ) + B 1 + --- , (153) 

with {X, Y} being locally defined Cartesian coordinates as in the a = 1/2 case. 
Again, if B 2 7^ 0, then on the 2-dimensional surface, E, does not correspond to a 
solution to the equation of interest but rather to an equation with an additional 
5-function source. Thus, again the solutions with B 2 7^ are not acceptable. 
Thus, hereafter we discard the second solution F(cos 2 :r) in all cases. 

Having set B 2 = 0, we now turn our attention to the boundary at x — 0, 
corresponding to infinity in AdS spacetime. In order to more easily see the 
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asymptotic behaviour of <£> near the boundary, it is convenient to transform the 
argument of the hypergeometric function F from cos 2 x to sin 2 x. Using the 
linear transformation formulae for hypergeometric functions [32], we obtain the 
following expressions: 
For v 7^ 0, v ^ N, we have 

$ = B l • (sina:)^ 1 / 2 • (cosx) CT+1 / 2 • F {(» a , 1 + *5 ^ *) 

= B, ■ G v {x) ■ { lf c l^l ■ F (a,., C^, 1 - v, sin 2 x) 

+ nc- + Jnc:!) ' {sinx)2u ' F (c - ? c -' 1 + " ; sin2x) } ' (154) 



where G„(x) was defined in eq. (130) above. 
For v E N, we have 

$ = 5 1 -(sinx)^ 1 / 2 -(cosxr +1 / 2 -F(C^,C7,l + ^;cos 2 x) 

r(i + a)r(i/) ^ (CV)*(C2> ^. 



= B 1 -G v (x) 



7.1/1 ,^ ( Silla; ) 

fc=0 

(-i)T(i + o-) 2 ^ (C^ ff )fc(C^)* /. n2/« f, /■ 2 \ 

— -G— - . • (sinx) • > — rb • smi • < logfsm x) 



r(C^)r(c-)^ fe!(i-i/) fc 

oo 



+ 1) - ik* + 1/ + 1) + v(c + fc ) + + fc )} 

For z/ = 0, we have 

$ = J B 1 -(sin a :) 1 / 2 .(cos a :) CT+1 / 2 -F(Co^,Co7 ) l + ^cos 2 a;) 
d /-» /■ \ r(l + cr) (Cv)fc(Co,ff)fc / . N2fc 



(155) 



X 



[ 2 ^{k + 1) - V(Co" CT + fc) - ^(Co7 + fc ) - log(sin 2 x)} . (156) 



When v ^ 1 (i.e., case (i) above), we see from eq. (154) for v £ N (and eq. (155) 
for i/GN) that $ fails to be square integrable, except for the case in which u is a 
real number such that either or becomes zero or a negative integer, i.e., 

uj = T (2m + l + is + a), m = 0, 1, 2, .... (157) 

(When eq. (157) holds, then either r(^ CT ) or T((~^) diverges, so the first term 
of the right side of eq. (154) for v ^ N (or the first term of eq. (155) for v G N) 
vanishes.) When < v < 1 (case (ii)) we see from eq. (154) that $ is square 
integrable for all u G C. When v — (case (iii)) we see from eq. (156) that $ 
also is square- integrable for all w G C. 
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3.3 Self-adjoint extensions 

In this subsection, we shall determine all of the positive self-adjoint extensions 
of the operator A, eq. (128) that satisfy the regularity condition at the origin, 
x = 7r/2, discussed in the previous subsection. (As discussed in the previous 
subsection, this additional regularity condition need be imposed only in the case 
of the I = modes of a scalar field in 3 and 4 spacetime dimensions.) We have 
already seen in Proposition 3.1 that if v 2 < (case (iv)), A itself is unbounded 
below, so, clearly, no positive self-adjoint extensions exist. Thus, we need only 
consider cases (i), (ii), and (iii). In these cases, A is positive, so there always 
exists at least one positive self-adjoint extension. 

As explained in subsection 3.1 and at the beginning of subsection 3.2, we 
will determine all of the possible self-adjoint extensions of A by finding all of 
the square- integrable solutions to eq. (143) with oo 2 = ±2i. We have just seen 
that in case (i) (i.e., v 2 ^ 1), there are no such solutions, so the deficiency 
subspaces are trivial. Hence, there is a unique self-adjoint extension of A, which 
is automatically positive. Thus, dynamics is well defined in this case, without the 
need to specify any additional boundary conditions at infinity. One may interpret 
this result as being a consequence of the fact that the effective potential term in A 
is sufficiently repulsive near infinity that the fields are completely reflected back 
without ever reaching the boundary. Hence, no boundary conditions at infinity 
need be imposed. 

On the other hand, in cases (ii) and (iii), we found in the previous subsection 
that there is precisely one linearly independent square integrable solution to the 
eigenvector equation (143) with u 2 = ±2i that satisfies the regularity condition 
at the origin. Consequently, each of the deficiency subspaces J£± (for the operator 
2A) is one-dimensional, and is spanned by the vector <&± of eqs. (154) (for case 
(ii)) and (156) (for case (iii)), with uj = l±i, respectively, and with (-£>i)± chosen 
so that ||$+||l2 = ||$-IU 2 = 1- 

Since n + = n_ = 1 the only (partial) isometries from JtT + to j£l are given by 

U a <Z> + = e ia $_ (158) 

where a G (— 7r, 7r]. Hence, in cases (ii) and (iii) there exists a one-parameter 
family of self-adjoint extensions, A a , of A, (see eqs. (139) and (140)) parametrised 
by a G (— 7r, 7r]. Our remaining tasks are (1) to interpret the choice of extension, 
A a , in terms of a choice of boundary condition at infinity and (2) to determine 
for what values of a the extension A a is positive. 

In order to investigate both of these issues, we need to know the asymptotic 
behaviour of 

$ a = $+ + e ia $_ (159) 
near x — 0. For < v < 1, we obtain 

$ a oc G„(x)-{a u + K(smx) 2u + c„(smx) 2 + ---} , (160) 
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where G v (x) was defined in eq. (130) above. Here the coefficients of the two 
leading terms are given by 



2e ia/2 T{l + a)F{ 



a-(I-H) 



K = 



r 

2e ia / 2 T(l + a)T(-is 



sm 



r ( cs: i) 



a 

tt • sm ' 
2 v 2 



(161) 
(162) 



where we have used the facts that = (, 
_ A -(i+0« 



-(1+0* 



+ i, c 



1-i 



+ 1, and 



, and where we have defined 6± v £ (— 7r/2,7r/2] by 
sm 0±„ = 



(163) 



y/1 + (±i/ + a) 2 ' 
Similarly for z/ = we obtain 

$ a oc G (x) • |a log(sin 2 x) + 6 + c (sinx) 2 log(sin 2 x) H | . (164) 

Using the property of the digamma function, ip(z+l) = ip(z) + l/ z (see eq. (149)), 
we find that the leading two coefficients are 



2e W2 r ( 1 + a ) 



a 



bo 



(l-H) 



(C +a ) 



'a 

sm i — - 0, 



(165) 



2e ia / 2 r(l + a) 



Co 



-(i+0 



r C 



(i+0 



+2 | 7 + Re [V (CoTi 



(i+0 



2 cos 2 # ' cos ^ — 

+ sin # cos #o | • sin — 9 ^j 



with 7 being Euler number and with 9 G [0, 7r/2) defined by 



smt/o = 



, (166) 



(167) 



By inspection of eqs. (161)-(162) we see that in case (ii) (i.e., < v < 1), the 
value of a is in one-to-one correspondence with the value of the ratio b u /a u , which 
may take any real value (including ±oo, corresponding to a v — 0). Thus, in case 
(ii), the choice of self-adjoint extension of A is uniquely specified by the value 
of b u /a u , and a self-adjoint extension exists for any specified real value of this 
quantity (including ±oo). Similarly, inspection of eqs. (165)-(166) shows that in 
case (iii) (i.e., v — 0), the self-adjoint extensions of A are uniquely characterised 
by the ratio bo/ao, which may assume any real value (including ±oo). 
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The relationship between the choice of extension and the choice of boundary 
conditions can now be seen as follows. It is easily seen that for arbitrary initial 
data in the general evolution law (4) yields <pt G Dom^e). However, from 
eq. (139), we see that Doul(Ae) consists of vectors in the domain of the closure 
of A — which have as rapid fall-off as is possible at the boundary — together with 
vectors of the form $ a = $ + + Therefore, the boundary conditions that 

are satisfied by vectors of this latter form are the boundary conditions that will 
be satisfied by all solutions that arise from arbitrary initial data in C^°. Con- 
sequently, the choice of self-adjoint extension in cases (ii) and (iii) corresponds 
precisely to the specification of the allowed asymptotic behaviour (as determined 
by the values of b v ja v and b /a , respectively) of solutions near infinity. 

The interpretation of the allowed boundary conditions at infinity is most ap- 
parent in the conformally invariant case, v = 1/2, since in this case, the effective 
potential term in the operator A is regular at the boundary. For v = 1/2, it is 
easily seen from eq. (160) that we have 



d§ a /dx 



a x=0 



h " (168) 



Therefore, we see that the choice b u /a u = ±oo (i.e., a v = 0) corresponds to the 
Dirichlet boundary condition, § a \x=o — 0, and the choice b v ja v = corresponds 
to the Neumann boundary condition, d$ a /dx\ x= Q = 0. The other real choices 
of b v ja v correspond to general Robin boundary conditions, in which a linear 
combination of 3> a | x =o and d§ a /dx\ x=0 is required to vanish. 

For v 7^ 1/2, the effective potential term in A is divergent at x = and the 
ratio (d& a /dx)/Q a itself is no longer well-defined at x — 0. Nevertheless, for < 
v < 1, the dominant behaviour of G~ x <& a as x — > is governed by a u , whereas the 
dominant behaviour of d(G~ 1 ^ a )/dx is governed by b v . Therefore, for < v < 1, 
it seems natural to define "generalized Dirichlet conditions" as corresponding 
to the case a u = and "generalized Neumann conditions" as corresponding to 
the case b v = 0. The other choices of self-adjoint extensions may be viewed 
as corresponding to generalized Robin conditions 5 . The Friedrich's extension 
corresponds to the choice a u = 0, i.e., to generalized Dirichlet conditions. 

For v — 0, the dominant behaviour of both G 1( l ) a and d(G^ 1 § a )/dx as 
x — > is governed by ao- Thus, we may interpret the self-adjoint extension with 
do = as corresponding to the simultaneous imposition of generalized Dirichlet 
and Neumann conditions. This extension is precisely the Friedrich's extension. 

5 Generalized Robin boundary conditions play a role in the AdS-CFT correspondence when 
one considers the double-trace perturbations of the quantum field theory on the boundary at 
infinity. In this context, the ratio b v ja v appears in the action for the boundary field theory 
as a coupling constant to the double-trace interaction, and generalized Dirichlet and Neumann 
conditions correspond to fixed points of a renormalization group flow in the quantum field 
theory on the boundary [12] (see also, e.g., [33, 34, 35, 15] and references therein). 
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We turn now to the determination of which of the self-adjoint extensions A a 
are positive. Our main results are given by the following theorem: 

Theorem 3.2: Let A be the operator eq. (128) on the Hilbert space L 2 ([0, 7r/2]), 

defined on the domain C^°(0,7r/2). Then the self-adjoint extensions of A — 

satisfying the regularity condition 6 at the origin x = 7r/2 (see subsection 3.2)— 

and the positivity properties of these extensions are as follows: 

case (i) If v 2 ^ 1, then there exists a unique self-adjoint extension, and it is 

positive. 

case (ii) If < v 2 < 1, then the self-adjoint extensions of A comprise a one- 
parameter family, A a . These extensions are characterised by the ratio b v ja v , 
which may assume any real value, but the extensions are positive if and only if 
b v ja v satisfies 



T(-v) 



r u) 



2 ' 



(169) 



(This inequality includes the case b v /a u = oo, i.e., generalized Dirichlet condi- 
tions, a v — 0.) In terms of a, this condition is expressed as 



sin (a/2 - 6_ 
sin (a/2 — 9 +v ) 



'! + (-// + of 
\ + {u + a) 2 





2 r(C»„) 2 


r (ti M> ) 


2 r(C) 2 



(170) 



case (iii) If v 2 = 0, then the self-adjoint extensions of A comprise a one- 
parameter family, A a . These extensions are characterised by the ratio &o/ a o> 
which may assume any real value, and are positive if and only if 



^ < 2 7 + 24 (C) • 
do 



(171) 



(This inequality includes the case b /a = — oo, i.e., generalized Dirichlet-Neumann 
conditions, ao — 0.) In terms of a, this condition is expressed as 



cos (a/2 — 9 ) 
sin (a/2 - O ) 



^(l + ( 7 2 )-U 



(7+1 



Re 



a — % 



a 



1 + a 2 



• (172) 



case (iv) If v 2 < 0, then A itself is unbounded below, so all self-adjoint extensions 
of A are unbounded below. 

Proof: We have already proven the above claims for cases (i) and (iv), and we 
have also established that in cases (ii) and (iii), the self-adjoint extensions of A 

6 As discussed in subsection 3.2, this regularity condition is relevant only for the / = modes 
of a scalar field in 3 and 4 spacetime dimensions. 
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comprise a one-parameter family, A a . Therefore, it remains only to prove the 
positivity claims in cases (ii) and (iii). 

It is worth noting that by a relatively straightforward calculation A a can be 
seen to be positive whenever b u /a u ^ in case (ii) and when ao = in case (iii). 
Namely, to prove the positivity of A a , it suffices to prove its positivity on the 
smaller domain V a defined by 

V a = {% + c$ a |$ G Dom(A),c G C} , (173) 

since it is clear that A a is the closure of its restriction to T> a and positivity is 
preserved under taking closures. Note that all $ G T> a are smooth. Furthermore, 
since A a (& + c$ Q ) = A$ + 2ic$ + — 2ice"*<E>_ (where the factor of 2 difference 
from eq. (140) arises from our present use of ±2i rather than ±i in the definition 
of the deficiency subspaces), it follows that for all $ G T> ai we have 



J \ ox 2 sin x cos 2 x 

As in Proposition 3.1, we may write 



A = D~D + + (-u + a + l) 2 (175) 

with 



where G v (i) is given by eq. (130). 

In case (ii) — i.e., when < v < 1 — it can be checked that _D + $ is square- 
integrable for all $ G V a . Consequently, integration by parts can be justified, 
and we obtain 

($,A$) L 2 = [-$* J D+$r /2 + (-v + a+ l) 2 / ' |$| 2 rfx+ f 1 \D + $\ 2 dx.(l77) 

The boundary term from the origin, x = tt/2, vanishes. The boundary term from 
x = is easily evaluated to yield 

($, A a $) L 2 = 2va*X + (-u + a + l) 2 ||$||i 2 + \\D + n 2 L 2 , (178) 

where a u and b u are given by eqs. (161) and (162). Note that va* v b v is always real 
(since b v /a u is real for any self-adjoint extension), as must be the case, since A a 
is self-adjoint. By inspection, A a will be positive whenever a* v b u ^ 0, as is the 
case whenever b u /a u ^ 0, including the case a v = 0. 

In case (iii) — i.e., when v = — we observe from (164) that the asymptotic 
behaviour of D + <& a near x = is 

D + $ a = 2(cosa;) a+3/2 {a (sina;)" 1/2 + c (sina;) 3/2 log(sin 2 :r) + •••}. (179) 
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We see, therefore, that D + Q a is not square integrable except for the case when 
a = 0. In this latter case, integration by parts is again justified, and we have 

($,A Q $) L2 = (<r+ l) 2 ||$||i 2 + \\D + m 2 > 0, (180) 

so A a is positive, as claimed. 

To determine the precise range of positivity of A a in cases (ii) and (iii), we 
note that since A a is strictly positive on the original domain of A (namely, the C£j° 
functions) and since the deficiency subspaces are only 1-dimensional, the negative 
spectral subspace of A a can be at most 1-dimensional. If the strictly negative 
spectral subspace of A a is 0-dimensional, then A a is positive. On the other hand, 
if the strictly negative spectral subspace of A a is 1-dimensional, then A a must 
have an eigenvector, typ, of negative eigenvalue, — (3 2 , with f3 > 0. By elliptic 
regularity, ^p must be a smooth solution to the differential equation (143), 

^ + ^ + ^W = -^. (181 ) 



dx 2 sin x cos 2 x 

Furthermore, in order to lie in the closure of T> a , it is necessary that the asymp- 
totic behaviour of ^ p approach that of c§ a at x = for some constant c G C in 
the sense that as x — > we must have ^ p — c& a = o([sm(x)] u+1 / 2 ) . 

In case (ii), it follows from our general solution, eq. (154), with u = i/3 that in 
order to have the required asymptotic behaviour of ^ p as x — > 0, (3 must satisfy 

ir(or Km (l82) 



r (CQ 



a u T(-u) ' 



It follows from the properties of the gamma-function that for fixed a ^ and fixed 
v in the range < v < 1, the left side of this equation can be made arbitrarily 
large by letting j3 — > oo and takes its minimum at f3 — 0. Taking into account 
the fact that T{v)/T{—v) is negative for < v < 1, it follows that a solution to 
eq. (181) with f3 > satisfying the required asymptotic behaviour exists if and 
only if 



r(-u) 



r(iz) 



r(c° ) 2 

j . (183) 



1 



Thus, eigenvector of negative eigenvalue exists — and hence A Q fails to be positive — 
if and only if this equation holds. This proves the claims of the theorem for 
case (ii). 

Similarly, in case (iii), it follows from our general solution, eq. (156), that in 
order to have the required asymptotic behaviour as x — > 0, (3 must satisfy 



2 7 + ^(c£)+^(C)=^" (184) 
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From the properties of the digamma function, ip(x), we see that the left side 
achieves its minimum value when /3 — and becomes arbitrarily large when 
/3 — > oo. Thus, we can solve this equation — and thereby obtain an eigenvector of 
A a of negative eigenvalue — if and only 

^>2 7 + 2^(Co%) • (185) 
a 

This proves the claims of the theorem for case (iii). □ 
We conclude with several remarks. 

As already noted, case (ii) was considered by Breitenlohner and Freedman [8, 
9] for scalar fields in 4-dimensional AdS spacetime. They showed that with the 
imposition of generalized Dirichlet or Neumann boundary conditions, the initial 
value problem is well-posed in this case. Their choice of these boundary conditions 
arose from the requirement that an energy functional for the scalar field be finite, 
positive, and conserved. However, they considered only two types of energy 
functionals: the usual energy functional of a minimally coupled scalar field, and 
a conformal ("improved") energy functional. The requirement that the usual 
energy functional be conserved allows only the generalized Dirichlet boundary 
condition, whereas conservation of the conformal energy allows us to take either 
the generalized Dirichlet or the generalized Neumann boundary condition. We 
wish to emphasise here that for any choice of a positive self-adjoint extension Ae, 
it is always possible to define a conserved, positive energy E for all $ G ~Dova(A E ) 
by 

E($) = + ($,A E $)v, (186) 

where $ = d<&/dt. The generalized Dirichlet and Neumann conditions considered 
in [8, 9], are merely two values of the one-parameter family of choices given by 
eq. (169). For each boundary condition in this range, there is a corresponding 
conserved, positive energy (see also [11] for another definition of such a positive, 
conserved energy in the scalar field case). 

As we have proven above, in cases (ii) and (iii), A also admits non-positive 
self-adjoint extensions; in case (iv), all self-adjoint extensions are unbounded 
below. If one chooses a non-positive self-adjoint extension, then one can still 
define dynamics in a manner analogous to eq. (4), but the dynamics will be 
unstable in the sense that generic solutions will grow unboundedly in time (see 
[36]). Indeed, in case (iv), since A is unbounded below, generic solutions will 
always grow in time more rapidly than exp(at) for any positive a. 

We note that our analysis shows that in AdS spacetime of any dimension, it 
is always possible to define a stable dynamics for electromagnetic and gravita- 
tional perturbations. In spacetime dimension ^ 7, the dynamics is unique, and no 
choice of boundary conditions at infinity need be made at all. However, in 4 and 6 
dimensions, one must make a choice of boundary conditions (in the range defined 
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by eq. (169)) for the scalar type of electromagnetic and gravitational perturba- 
tions; in 4 dimensions, one must also choose boundary conditions in this range for 
the vector type of electromagnetic and gravitational perturbations. Similarly, in 
5-dimensional AdS spacetime one must make a choice of boundary conditions (in 
the range defined by eq. (171)) for the scalar type of electromagnetic and gravi- 
tational perturbations. Each of these choices of boundary conditions corresponds 
to a different theory of classical bulk dynamics. 

Particularly interesting in the context of the AdS-CFT correspondence is the 
case in which the bulk geometry contains a black hole [37]. In the Schwarzschild- 
AdS background, we also can obtain master equations for perturbations [30] 
analogous to those derived in this paper. We then expect the same classification 
(i)-(iv) of the corresponding operator A to hold, because the classification (i)-(iv) 
is made by inspecting merely the asymptotic behaviour of perturbations near the 
boundary. However, since the effective potential term in the corresponding oper- 
ator A will change significantly near the black hole event horizon [38], we would 
expect that the precise relation between the positivity of self- adjoint extensions 
and the boundary conditions at infinity will change. 
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Appendix 

We show here that for the / = 1 modes, there exist no dynamical degrees of 
freedom for gravitational perturbations (i.e., no dipole gravitational radiation). 
(See also [28].) 

Consider, first, the vector type of perturbations. For I = 1 (k v = n — 1), V, 
becomes a Killing vector field on the n-sphere. Then 2D^h^^ in (73) vanishes, 

and therefore is not defined. In this case, we can find a gauge-invariant 
variable 
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Z (i) 



r 2 V, 



( 



H (i) 




-r 2 V b 



( 



H (i) 




(187) 



ab 



a 
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which has only a single independent component. Then the linearised Einstein 
equations for Z^ ab become 

V a (r n Z^ ab )=0, (188) 
and can be explicitly solved, yielding 

Z (1) ab = e ab ^ , (189) 

where C is a constant. Therefore the vector- type perturbation of the I = 1 
mode is a non-dynamical (i.e., time independent) perturbation, which can be 
interpreted as describing a change in the angular momentum of the background. 

For I = 1 (kg = n) of scalar-type perturbation, £)j§ is a conformal Killing 
vector field, that is, DiDjS* + 7jjS identically vanishes. Hence the term hj) in 
(73) is not present. Setting = in the definitions (99) and (100) (and 
X a = —H(°\), we can still use Z and Z ab as metric perturbation variables, but 
in this case, Z and Z ab are no longer gauge-invariant and equation (103) is not 
obtained. Under gauge-transformations generated by ^°\E>dy a + DiExiz 1 , we 
find in particular, Z a a - Z -> Z a a - Z + 5(Z a a - Z) with 

- Z) , \V<% - <"- 2 »'"- 4 » l - j E , (190) 



r 



where S = — 2r^ n_4 ^ 2 ^ ^. This implies that we can make eq. (103) hold by choos- 
ing S suitably, or in other words we can impose eq. (103) as a gauge condition. 
In this way, even for the / = 1 modes, we can recover the same expression of the 
linearised Einstein equations for Z and Z ab as in the / ^ 2 case. 

We should note however that the gauge condition Z a a — Z = does not com- 
pletely fix gauges; it is possible to take further gauge transformations generated 
by S that satisfies the homogeneous equation 

^_("-2)("-4)l_n(n + 2)l| 5 = ^ (m) 



r 



f: 2 



Therefore the remaining gauge degrees of freedom are characterised by two arbi- 
trary functions of t, i.e., initial data (E,dE/dr) on a Cauchy surface r = r in 
the globally hyperbolic spacetime (O 2 , —g a b)- 

Now we shall show that the solutions Z and Z ab to the linearised Einstein 
equations are gauge modes; Z and Z ab can be made to vanish by further gauge 
transformations. First, with the gauge condition (103), repeating the same ar- 
gument as in the / ^ 2 case (below eq. (104)), we can express Z ab in terms of a 
scalar function as 



V V 6 - ^<7o6 ) o, . (192) 
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where 0^ satisfies the linearised Einstein equation (see eqs. (120) and (121) with 
k 2 s = n) 



V a r ~ n — 2 

£(<^) = r A [ V a V a - n— -V a + ) & = V(c 2 sin t + c 3 cost) , (193) 



with C2, C3 being arbitrary constants. Note that the term c'/r, which may exist in 
the right hand side (see eq. (123)), has been already set to zero by using the resid- 
ual freedom of 0^ defined by eq. (118). Next, we find the gauge transformation 
z ab -> Z' ab = Z ab + 5Z ab with 



5Z ab = 



„ (n-2) V a r~ „ V b r~ „ 



-V 6 ~ + — V„~ 
2 \ r r 



, . V c r e „ n(n-2)(V a r)V fe r„ 
+ (n - l)# a & V c ^ + 



nV a V b r„ 



+(n - l)g, 



ab 



1 n (V c r) V c r 



(194) 



where S satisfies eq. (191). To achieve the gauge = 0, we consider on an 
initial Cauchy surface r = r the equations 



(195) 



^ 6 = 5Z ab (~) + ( V a V b - -g ab ) & = . 



The trace part of eq. (195) is expressed as 

d ( 1 



y r (n-2)/2 



1 



V d fr n/2 



nr 



(n+2)/2y 



f n /2 Q r \ y 
(c 2 sin t + c 3 cos t) , 



(196) 



where = r~ n / 2 0£ and we have used eqs. (191) and (193). Note that the 
constants c 2 and C3 in the right side of this equation are actually required to 
vanish due to the consistency with the (t, r)-component of eq. (195), and that 
in the case obeys the homogeneous master equation E(r n//2 $ ? ) = 0. The 
(t, t) -component of eq. (195) is expressed as 



ldV 2 . ^V 2 



dr^ 



V^dt 2 ^ 



n dV 2 (ra-1) n(n-l)V 2 \^ 



4r dr r 2 
1 d 2 1 dV 2 d 
V^dt 2 + 2~dr~dr~ + 



+ 



2 r 
w <9\/ 2 _ 1_ 

47^7 ~~ Z 2 



* ■ 



(197) 
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For any (<3>£, d&^/dr), solving eqs. (196) and (197), we obtain the initial data 
(S, dS/dr) that satisfy eq. (195) on the initial surface r = r . Then we define S 
to be the solution of eq. (191) with these initial data so that Z' ab = is achieved 
throughout (0 2 ,-g ab ). 
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